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This  study  concerns  the  mathematical  analysis  of  constrained-layer  damping  of  extensional 

waves  in  plates  of  infinite  extent,  with  and  without  fluid  loading.  Previous  work  was  _ 

mostly  limited  to  flexural  waves.  Some  aspects  of  fluid  loading  for  flexural  waves  may  be  ;ributlon( 

understood  by  means  of  thin-plate  theory.  Therefore,  a  similar  theory  was  developed  for  - 

extensional  waves.  The  description  and  examples  presented  here  are  based  on  three  Availability  Codes 

models:  The  first  is  an  extension  of  Kerwin’s  1959  model  [E.  M.  Kerwin,  J.  Acoust.  Soc. 


Am.  31,  952-962  ( 1959)],  the  second  is  a  hybrid  model  in  which  the  base  plate  is  described  by  t 
exact  elasticity  theory  and  the  other  two  layers  by  Kerwin’s  concepts,  and  the  third  uses 
exact  elasticity  theory  for  all  three  layers.  It  is  shown  that  the  extended  Kerwin  model  is  useful 
in  the  design  of  constrained-layer  damping  for  extensional  waves  as  well  as  for  flexural  \ 
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INTRODUCTION 

The  technique  of  constrained-layer  damping  (CLD) 
of  acoustic  waves  consists  of  attaching  a  thin  elastomer 
layer  with  high  viscoelastic  loss  plus  a  stiff  covering  layer 
to  a  bar,  plate,  or  structure.  This  stiff  layer  forces  the  elas¬ 
tomer  into  shear,  with  concomitant  large  loss,  as  compared 
with  purely  extensional  damping  loss  in  the  elastomer  layer 
without  the  constraining  layer.  Although,  for  a  soft  elas¬ 
tomer,  the  loss  tangent  of  the  shear  modulus  G  is  about  the 
same  as  the  loss  tangent  of  the  Young’s  modulus  E,  the 
energy  loss  in  the  constrained  layer  is  of  the  order  l/khi  of 
the  loss  in  a  comparable  unconstrained  layer,  where  k  is 
the  wave  number  in  the  composite  plate  and  /12  is  the  thick¬ 
ness  of  the  elastomer  layer. 

This  physical  explanation  of  the  effect  is  represented  in 
a  model  described  by  E.  M.  Kerwin  in  a  classical  paper. '  In 
this  model  the  base  plate  is  treated  by  thin-plate  theory, 
thus  for  small  values  oi  khi,  where  Aj  is  the  thickness  of  the 
base  plate.  His  model  gives  good  results  within  the  appli¬ 
cable  restrictions;  flexural  waves  at  low  frequency  in  a  base 
plate  or  bar  with  thin  added  layers.  Extensional  waves  are 
not  considered,  and  fluid  loading  is  not  accounted  for. 

It  is  typical  fer  Kerwin’s  model  that  the  attenuation 
constant  a  is  given  as  an  explicit  algebraic  expression  in 
terms  of  the  physical  and  geometric  parameters  of  the  com¬ 
posite  structure.  This  is  an  advantage  in  studying  the  in¬ 
fluence  of  the  various  parameters  on  the  attenuation.  This 
author^  attempted  to  construct  an  extension  of  Kerwin’s 
model  that  retained  this  advantage,  white  relievmg  some  ot 
the  restrictions.  A  short  description  of  the  extended  Ker¬ 
win  model  is  given  in  Appendix  A,  to  serve  as  a  basis  for 
the  examples  in  the  following  text.  In  order  to  include 
extensional  waves,  it  was  necessary  to  account  for  inertia  of 
the  constraining  layer.  Moreover,  it  became  clear  that  the 
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results  of  the  extended  Kerwin  model  were  greatly  im¬ 
proved,  when  wave  speeds  derived  from  thick-plate  theory 
were  used  instead  of  the  thin-plate  speeds  used  by  Kerwin. 
This  was  true  for  both  flexural  and  extensional  waves. 

In  order  to  evaluate  the  properties  of  this  extended 
Kerwin  model  and  to  check  the  range  of  its  validity,  a 
model  following  from  exact  elasticity  theory  for  all  three 
layers  was  developed.  As  an  intermediate  step,  a  hybrid 
model  was  constructed  that  contains  features  of  both  the 
extended  Kerwin  and  exact  models.^  This  model  entails 
fewer  mathematical  operations  than  the  exact  model,  and, 
as  a  consequence,  it  will  produce  results  faster. 

The  two  types  of  freely  propagating  waves  in  an  infi¬ 
nite  plate  indicated  by  the  names  flexural  (antisymmetric) 
and  extensional  (symmetric),  collectively  named  Lamb 
waves,  are  mathematically  distinguished  by  the  parity  of 
the  functions  representing  the  field  variables  with  respect 
to  the  z  coordinate,  perpendicular  to  the  faces  of  the  plate, 
with  the  origin  halfway  between  the  faces.  It  follows  that 
these  two  types  can  each  occur  in  a  pure  form,  without 
admixture  of  the  other  one,  only  in  a  situation  that  is 
strictly  symmetric  with  respect  to  the  z  coordinate.  Thus 
the  plate  should  be  single  and  homogeneous  and  it  should 
be  symmetrically  loaded  by  the  same  type  of  fluid  on  both 
sides  (including  the  unloaded  case).  This  would  appear  to 
exclude  the  possibility  of  considering  either  one  of  these 
types  separately  when  the  plate  is  treated  by  constrained- 
layer  damping  on  one  side  only  and/or  under  one-sided 
fluid  damping. 

In  stuaying  the  results  from  exact  elasticity  theory  for 
fluid-loaded  plates  with  constrained  layer  damping,  when 
both  types  of  waves  are  accounted  for,  it  is  found,  though, 
that  these  two  types  of  waves  preserve  their  own  identity 
even  in  a  nonsymmetric  situation.  In  the  hybrid  and  fully 
exact  models  the  two  types  are  always  present  and,  at  low 
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values  of  kd,  one  finds  two  complex  roots  that  are  easily 
recognized  as  to  their  designation,  either  mostly  antisym¬ 
metric  or  mostly  symmetric.  In  computational  practice  one 
finds  one  or  the  other  wave  by  starting  the  iterative  root 
finder  from  a  seed  representing  either  type  of  wave,  com¬ 
puted  from  thick-plate  or  exact  theory  for  a  single,  not 
fluid-loaded  plate.  The  resulting  wave  speed  computed 
from  the  complex  root  is  affected  by  the  presence  of  the 
other  type,  but  not  to  a  large  extent  at  low  frequencies.  For 
high  frequencies  it  is  often  difficult  to  find  an  adequate 
seed.  In  experimental  realizations  one  may  create  the  de¬ 
sired  type  of  wave  by  the  method  of  excitation,  for  instance 
by  means  of  two  actuators  on  either  side  of  the  plate  driven 
with  opposite  or  equal  phase.  At  certain  values  of  the  fre¬ 
quency,  higher-order  modes  appear  on  the  scene,  deter¬ 
mined  by  thickness-shear  modes  for  flexural  waves  and 
thickness-dilatational  modes  for  extensional  waves,  with 
some  influence  by  the  other  type.  At  higher  values  of  kd, 
the  two  types  may  be  difficult  to  distinguish  since  the  roots 
are  close  together  and  corisiderable  accuracy  was  needed  to 
separate  them. 

The  extended  Kerwin  model  is  limited  by  the  fact  that 
a  wave  speed  for  a  purely  flexural  wave  or  a  purely  exten¬ 
sional  wave  has  to  be  entered  into  the  computation,  with¬ 
out  accounting  for  the  possible  effect  of  the  other  type.  The 
hybrid  model  uses  exact  elasticity  theory  for  the  base  plate 
and,  therefore,  the  wave  speed  and  attenuation  are  both 
part  of  the  root  finding  effort,  with  proper  account  of  the 
other  type,  in  asymmetric  situations.  It  shares  this  aspect 
with  the  fully  exact  model.  The  fully  exact  model  allows 
waves  with  such  a  large  value  of  kd  that  even  the  elastomer 
layer  and  the  constraining  layer  do  not  behave  like  “thin” 
layers.  Some  examples  of  this  are  given  in  this  paper;  the 
number  of  possible  waves  increases  greatly  and  exceeds  the 
range  for  which  the  notions  of  constrained-layer  damping 
are  practical.  Except  for  these  large  values  of  kd,  the  wave 
speeds  as  computed  by  the  fully  exact  model  do  not  differ 
significantly  from  the  thick-plate  or  exact  values  computed 
for  a  single,  unloaded  plate. 

The  following  outline  lists  the  three  models  for 
constrained-layer  damping  used  in  this  study,  with  their 
characteristic  properties. 

1.  Extended  Kerwin  model: 

(a)  Flexural  or  extensional  waves. 

(b)  Inertia  of  constraining  layer. 

(c)  Complex  shear  parameter. 

(d)  Wave  speed  derived  from  thick-plate  theory  for 
single  unloaded  plate. 

2.  Hybrid  model: 

(a)  Exact  elasticity  theory  for  base  plate. 

(b)  Other  layers  as  in  extended  Kerwin  model. 

(c)  Simultaneous  flexural  and  extensional  waves, 

(d)  With  or  without  fluid  loading. 

(e)  Wave  speed  computed  as  part  of  the  complex 
rootfinding  process. 

3.  Exact  model: 

(a)  Exact  elasticity  theory  for  all  three  layers. 

(b)  Simultaneous  flexural  and  extensional  waves. 


(c)  With  or  without  fluid  loading. 

(d)  Wave  speed  computed  with  proper  accouni  of 
nonsymmetric  fluid-loading  and  influence  of 
added  layers. 

In  a  previous  study,  the  predictive  capability  of  the 
extended  Kerwin  model  was  studied  for  the  case  of  flexural 
waves  and  the  effect  of  fluid  loading.^  Preliminary  results 
of  damping  of  extensional  waves  by  CLD  were  given  in 
short  form  in  Ref.  4.  The  present  paper  gives  an  expanded 
account  of  the  application  of  the  three  models  to  damping 
of  extensional  waves  in  infinite  plates  without  and  with 
fluid  loading. 

For  flexural  waves,  valuable  insight  into  the  effect  of 
fluid  loading  may  be  gained  from  thin-plate  theory  with 
fluid  loading.*  Therefore,  i  model  was  developed  and  pre¬ 
sented  in  this  study  for  extensional  waves  in  a  thin  plate 
under  fluid  loading.  Some  results  are  given,  based  on  a 
derivation  (Appendix  B)  of  this  model.  After  this  discus¬ 
sion,  this  paper  continues  with  predictions  from  the  three 
models  for  extensional  waves  in  constrained  plates  without 
and  with  fluid  loading.  It  is  shown  that  the  behavior  of  the 
attenuation  as  a  function  of  frequency  is  strongly  affected 
by  the  ratio  of  the  extensional  wave  speed  in  the  constrain¬ 
ing  layer  relative  to  that  in  the  base  plate. 


I.  THIN-PLATE  THEORY  FOR  EXTENSIONAL  WAVES 
IN  A  FLUID-LOADED  PLATE 


In  this  paper  it  is  assumed  that  a  straight-crested,  har¬ 
monic  wave  propagates  in  the  x  direction  parallel  to  the 
faces  of  the  plate,  with  time  and  space  dependence  accord¬ 
ing  to  exp  i{kx—o>t).  The  complex  wave  number  k  is  writ¬ 
ten  as  k=k'  —  ia,  where  a  is  the  attenuation  constant. 

In  Appendix  B,  the  following  dispersion  relation  is 
derived  for  extensional  waves  in  a  fluid-loaded  plate,  for 
thin-plate  theory.  This  equation  is  given  in  terms  of  a  com¬ 
plex  variable  r,  which  is  related  to  the  wave  number  k  in 
the  plate  by  {kQT)^  =  k^  —  kl,  where  kg  is  the  wave  number 
in  the  medium,  with  wave  speed  Cq: 

r^_nr^  +  T[l-(co/cp^]-n[l-(co/c^)']=0.  (1) 

Here,  is  the  thin-plate  extensional  wave  speed  in  the 
plate,  and  Cj  the  dilatational  wave  speed  in  the  plate  ma¬ 
terial.  A  characteristic  frequency  co,.  and  a  corresponding 
dimensionless  frequency  ft  are  introduced  by 
pQ)cy(C(/i)  and  p,  and  pg  are  the  densities  of 

the  plate  material  and  medium,  respectively;  and  d  is  the 
half-thickness  of  the  plate. 

The  characteristic  frequency  o)^  for  brass  in  water  is 
1.63 X  10*’  rad/s,  and  thus  for  the  (low-)  frequency  region 
where  thin-plate  theory  may  be  applied,  the  nondimen- 
sioiial  frequency  ft  is  very  small.  As  a  consequence,  tncie 
exists  a  small  real  root,  for  small  ft,  analogous  to  the  real 
root  for  the  corresponding  equation  for  flexural  waves,* 
namely. 


~  l-(co/fp)^ 


(2) 
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TABLE  I.  Material  parameten>  and  derived  quantities. 


A.  Brass 

B  Aluminum 

C  Silver 

Young's  modulus 

E 

104  GPa 

71  GPa 

78  GPa 

Poisson's  ratio 

V 

0.37 

0,33 

0,37 

Density 

P 

8500  kg/m‘ 

2700  kg/m' 

10  500  kg/m' 

Shear  wave  speed 

2113  m/s 

3149  m/s 

1647  m/s 

Extensional  wave 
speed  (thin) 

3765  m/s 

5462  m/s 

2933  m/s 

Characteristic 
frequency,  for 
10-cm  plate 

1.63  Mrad/s 

1.09  Mrad/s 

1,22  .Mrad/s 

Dim.  less  freq. 
k^  —  lttf/c^ 

d=0  05  m;  /=  10  kHa 

1.49 

1,00 

1.07 

FIG.  1.  Attenuation  factor  a  relative  to  the  real  part  of  the  wave  number 
k'  as  a  function  of  frequency  /,  for  a  single  10-cm-thick  brass  plate, 
loaded  by  water.  Solid  curves  are  computed  according  to  exact  elasticity 
theory,  dashed  curves  according  to  thin  plate  theory.  Curves  1  are  for 
water  on  both  sides,  curves  2  for  water  on  one  side. 


By  division  one  determines  the  quadratic  equation  for 
the  two  remaining  complex  conjugate  roots, 

T^  +  ar  +  b—O,  (3) 

where 


One  of  these  roots  leads  to  a  wave  number  with  positive 
attenuation  constant  a.  (The  other  root,  with  negative  at¬ 
tenuation  constant,  has  no  physical  significance.) 

For  small  values  of  ft  the  phase  speed  is  very  close  to 
Cp,  and  the  attenuation  a  relative  to  the  real  part  of  the 
wave  number  k'  is  approximated  by 


a  Poojd  (Cf/Cp) ^  -  (co/cj) ^ 
p,  Co  [(l-(co/c^)2]'''^ 


Viscoelastic  layer  (hypothetical) 

Real  part  of  shear  modulus  Gi  =  10  MPa 

Loss  tangent  of  shear  modulus  P  =  1 

Bulk  modulus  ^>=1  GPa 

Fluid,  water 

Density  Pi)=998  kg/m' 

Wave  speed  c,,- m/s 


It  app)ears  that  for  extensional  waves  the  wave  energy 
moves  to  the  side  of  the  fluid,  opposite  to  the  behavior  for 
flexural  waves.  ^ 

There  is  a  characteristic  difference  between  the  influ¬ 
ence  of  fluid  loading  in  flexural  versus  extensional  waves. 
The  phase  speed  of  flexural  waves  increases  with  frequency 
from  zero  to  the  Rayleigh  wave  speed.  Thus,  for  most 
cases,  there  will  be  a  frequency,  named  the  coincidence 
frequency,^  for  which  the  flexural  wave  speed  is  equal  to 
the  wave  speed  in  the  fluid.  This  plays  an  important  role  in 
the  occurrence  of  radiation.  In  contrast,  for  most  combi¬ 
nations  of  plate  material  and  fluid,  the  extensional  wave 
speed  will  be  larger  than  the  wave  speed  in  the  medium  for 
all  frequencies. 


In  Fig.  1  a  comparison  is  shown  of  the  relative  atten¬ 
uation  a/k'  for  a  single  fluid-loaded  plate,  according  to 
thin-plate  theory  (dashed  curves)  and  to  exact  elasticity 
theory  (solid  curves).  The  plate  is  10-cm-thick  brass.  The 
curves  indicated  by  1  correspond  to  two-sided  loading  by 
water.  One  may  observe  that  the  agreement  between  the 
two  models  is  excellent,  except  at  higher  frequencies.  Al¬ 
though  the  boundary  conditions  for  extensional  waves  can¬ 
not  be  satisfied  for  one-sided  fluid  loading  without  addi¬ 
tional  flexural  waves,  one  might  adapt  the  thin-plate  theory 
by  taking  half  the  value  of  the  total  fluid  stress  term  tf 
(equal  to  the  sum  of  the  stresses  on  the  two  faces  of  the 
plate),  as  given  for  two-sided  loading,  and  carrying  this 
value  for  <f  through  the  analysis.  The  results  are  shown  in 
Fig.  1  by  the  curves  marked  2.  It  shows  that  the  attenua¬ 
tion  for  one-sided  fluid-loading  thus  computed  does  not 
compare  well  with  the  exact-elasticity  result. 

One  may  observe  in  Fig.  1  that  the  two  solid  curves 
marked  1  and  2  differ  by  a  factor  of  two  over  most  of  the 
frequency  range.  At  the  high-frequency  end  the  two  curves 
merge,  thus  indicating  that  the  radiation  in  the  plate 
loaded  on  one  side  is  the  same  as  for  loading  on  both  sides. 


11.  DAMPING  OF  EXTENSIONAL  WAVES 

In  order  to  show  the  relative  merits  of  the  three  mod¬ 
els  described  before,  various  combinations  of  material  and 
geometric  parameters  are  chosen.  The  generality  of  the 
results  would  be  more  clearly  revealed  by  expressing  all 
parameters  in  dimensionless  form,  including  the  abscissa  of 
the  curves  shown.  The  representation  of  parameters  and 
results  in  dimensional  form,  though,  offers  a  more  concrete 
feeling  for  the  phenomena.  It  is  remarkable  that  dimen¬ 
sionless  representation  of  wave  speeds  or  wave  numbers  for 
Lamb  waves  in  unloaded  plates  is  very  similar  for  different 
materials;  the  graphs  depend  only  slightly  on  Poisson’s 
ratio  (see  Ref  6).  A  useful  expression  for  the  frequency  in 
dimensionless  form  is  the  product  k^,  where  k^  is  the  wave 
number  for  shear  waves  in  the  plate  material  and  d  the 
half-thickness  of  the  plate.  To  appreciate  its  relationship  to 
the  frequency,  the  value  of  k^  for  the  various  materials 
and  a  plate  of  0.1-m  thickness  is  shown  in  Table  I.  The 
wave  speeds  in  the  hybrid  and  fully  exact  models  are  in 
each  case  computed  from  the  complex  wave  number  found 
by  the  complex  root  finder.  The  wave  speeds  are  only 
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slightly  influenced  by  the  added  layers  and  by  fluid  loading 
at  small  values  of  k^.  Only  at  large  values  of  kd  do  sizable 
deviations  occur  (as  one  seen  in  Fig.  8).  At  these  large 
frequencies  the  situation  becomes  complicated  and  is  de¬ 
pendent  on  the  given  combination  of  materials  and  thick¬ 
nesses. 

The  study  of  the  algebraic  expression  for  the  extended 
Kerwin  model  shows  some  peculiar  structures,  as  is  dis¬ 
cussed  below.  These  same  features  show  up  in  the  hybrid 
and  fully  exact  models.  In  addition,  these  two  latter  models 
show  features  at  high  frequency  that  are  not  contained  in 
the  '"Ttended  Kerwin  model.  Insofar  as  a  physical  interpre¬ 
tation  is  found  (e.g.  concerning  Fig.  8),  these  phenomena 
are  considered  to  be  general  in  nature,  and  not  a  product  of 
the  chosen  combination  of  specific  parameters. 

A  major  point  in  the  development  of  the  extended 
Kerwin  model  is  the  introduction  of  inertia  of  the  con¬ 
straining  layer.  This  corresponds  to  the  appearance  of  the 
factor  [1  — (c/c3}^]  in  the  right-hand  side  of  Eq.  (Al)  of 
Appendix  A,  where  C3  is  the  extensional  wave  speed  in  the 
constraining  layer,  at  low  frequency  given  by  the  thin-plate 
extensional  speed,  and  c  is  the  pertinent  wave  speed  in  the 
base  plate.  For  flexural  waves  at  low  frequency  it  is  gener¬ 
ally  true  that  c4,c-^,  and  this  factor  is  then  close  to  one.  For 
extensional  waves  at  low  frequency,  c  is  equal  to  Cp,  the 
thin-plate  extensional  wave  speed  of  the  base  plate  mate¬ 
rial.  Thus  it  follows  that,  at  low  frequency,  extensional 
waves  are  not  damped  by  the  constrained-layer  technique  if 
the  constraining  plate  is  made  of  the  same  material  as  the 
base  plate,  since  then  this  factor  reduces  to  zero.  Further, 
the  maximum  value  of  the  third  factor  in  the  right-hand 
side  of  Eq.  (Al),  given  by  Eq.  (A3),  is  dependent  on  the 
algebraic  sign  of  [1  — (c/C3)^].  It  will  be  shown  that  this 
leads  to  a  considerable  increase  of  the  maximum  damping 
for  extensional  waves,  if  the  extensional  wave  speed  in  the 
constraining  layer  is  less  than  that  in  the  base  plate,  other 
factors  being  equal. 

In  the  first  set  of  examples,  the  materials  and  layer 
thicknesses  are  identical  to  those  in  the  previous  work  on 
flexural  waves. There  the  base  plate  was  chosen  to  be 
brass,  the  elastomer  layer  was  given  hypothetical  values, 
independent  of  frequency,  and  the  constraining  layer  was 
aluminum.  The  physical  parameters  of  these  and  other  ma¬ 
terials  used  in  the  followiag  examples  are  given  in  Table  1. 
The  thickness  of  the  base  {^ate  is  10  cm  in  all  cases.  Equa¬ 
tion  (A2)  of  Appendix  A  shows  that  the  maximum  damp¬ 
ing  is  found  for  a  given  value  of  the  real  part  of  the  shear 
parameter  g',  depending  on  the  wave  speeds  in  base  plate 
and  constraining  layer  and  the  loss  tangent  of  the  shear 
modulus  of  the  elastomer.  According  to  the  definition  of  g 
this  maximum  can  be  moved  to  a  desired  frequency  by  the 
choice  of  the  product  of  the  thicknesses  of  elastomer  and 
constraining  layers  A2A3.  The  second  factor  in  the  right- 
hand  side  of  Eq.  (Al)  shows  that  the  damping  increases 
with  increasing  thickness  ^3.  Of  course,  there  are  practical 
limitations  to  a  large  hj,  apart  from  the  restriction  to  a  thin 
constraining  layer,  implicit  in  the  model.  The  ratio  of  the 
thickness  of  the  constraining  layer  hy  to  that  of  the  elas¬ 
tomer  layer  is  set  equal  to  2  (arbitrarily).  The  parame- 
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FIG.  2.  Attenuation  factor  a  relative  to  the  real  pan  of  the  wave  number 
A:'  as  a  function  of  frequency  /,  for  a  constrained  plate  in  vacuum.  The 
physical  parameters  for  the  three  layers  are  given  in  Table  I.  The  base 
plate  is  10-cm  brass,  the  elastomer  layer  is  1.24  mm  thick,  the  constrain¬ 
ing  layer  is  2.48-mm  aluminum.  The  solid  curve  is  computed  according  to 
the  exact  model,  the  dashed  curve  according  to  the  hybrid  model,  and  the 
dotted-dashed  curve  according  to  the  extended  Kerwin  model. 

ters  were  chosen  such  that  the  maximum  damping  was 
near  1000  Hz,  for  flexural  waves.  Then  ^2=  l  ^'^  mm,  and 
A3 =2.48  mm.  This  set  of  geometric  parameters  is  also  used 
for  the  first  set  of  extensional  wave  examples,  denoted  as 
case  A.  To  demonstrate  the  application  of  the  extended 
Kerwin  model,  the  next  case,  B,  employs  thicknesses  of  the 
constraining  and  elastomer  layers  that  produce  maximum 
damping  of  extensional  waves  near  1000  Hz.  Again  de¬ 
manding  a  ratio  of  2,  one  finds  that  A2=7.75  mm  and 
A3=15.5  mm. 

til.  MAXIMUM  DAMPING  AT  6200  HZ  (CASE  A) 

A.  Constrained  plate  In  vacuum,  aluminum 
constraining  layer 

In  Fig.  2,  a  threefold  comparison  is  shown  for  the 
relative  attenuation  constant  a/k'  for  extensional  waves 
propagating  in  a  constrained  plate,  without  fluid  loading, 
for  the  three  models.  The  physical  parameters  for  the  three 
layers  are  given  in  Table  I.  They  are  the  same  as  those  for 
the  flexural  wave  examples  in  Ref.  3.  The  thicknesses  of  the 
elastomer  and  constraining  layers  were  chosen  such  that 
the  maximum  relative  attenuation  for  flexural  waves  was 
near  1000  Hz,  while  hy=2h2- 

One  sees  that  the  maximum  relative  attenuation  for 
extensional  waves  has  shifted  to  about  6200  Hz.  This  is 
due,  first,  to  the  fact  that  the  numerator  in  the  right-hand 
side  of  Eq.  ( A2)  is  smaller  for  extensional  waves  than  for 
flexural  waves.  Second,  according  to  the  definition  for  the 
shear  parameter  g  in  Appendix  A,  it  is  inversely  propor¬ 
tional  to  the  square  of  the  wave  number  k'.  In  order  to 
maintain  a  given  value  for  g  at  the  larger  extensional  wave 
speed,  as  compared  with  the  flexural  wave  speed,  the  fre¬ 
quency  will  have  to  be  accordingly  higher. 

The  sharp  dip  in  attenuation  near  1 5  kHz  is  due  to  an 
“equivoluminal  mode,”’  whereby  the  tangential  velocity 
component  at  the  faces  of  the  plate  is  zero;  thus  no  shear 
exists  in  the  elastomer  layer,  and  thus  no  viscoelastic 
damping.  Of  course,  this  feature  does  not  show  up  in  the 
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FIG.  3.  Attenuation  factor  a  relative  to  the  real  part  of  the  wave  number 
k'  as  ^  function  of  frequency  /,  for  a  constrained  plate  loaded  by  water, 
according  to  exact  elasticity  theory.  The  physical  parameters  for  the  three 
layers  are  given  in  Table  I.  The  base  plate  is  10-cm  brass,  the  elastomer 
layer  is  1.24  mm  thick,  the  constraining  layer  is  2.48-mm  aluminum.  Solid 
curve:  water  on  both  sides;  dashed  curve:  water  on  constrained  side; 
dotted-dashed  curve;  water  on  opposite  side. 

extended  Kerwin  model,  which  does  not  account  for  any 
structure  of  the  wave  fields  in  the  base  plate.  Apart  from 
this  phenomenon,  there  appears  to  be  reasonable  agree¬ 
ment  between  the  results  of  the  three  models  in  the  middle- 
frequency  range. 

B.  Fluid-loaded  constrained  plate,  aluminum 
constraining  layer 

In  Fig.  3,  a  comparison  is  shown  for  the  relative  at¬ 
tenuation  constant  a/k’  as  a  function  of  frequency  for  a 
constrained  plate  loaded  by  water  on  both  sides,  on  the  side 
of  the  constraining  layer,  and  on  the  opposite  side,  com¬ 
puted  by  exact  elasticity  theory.  One  sees  that,  for  a  large 
frequency  range,  the  curves  are  practically  identical  to 
those  of  Fig.  1.  The  viscoelastic  damping,  which  is  much 
smaller  than  the  radiation  damping,  is  apparently  just 
added  to  the  latter,  without  much  influence  of  the  fluid  on 
the  composite  plate.  Only  at  the  high-frequency  end  is 
there  a  structure  different  from  the  one  in  Fig.  1,  which  is 
apparently  due  to  this  interaction:  the  two  curves  for  load¬ 
ing  on  both  sides  and  on  the  constrained  side  merge  and 
show  a  decrease  in  attenuation  not  observed  in  Fig.  1. 

C.  Constrained  plate  in  vacuum,  silver  constraining 
layer 

In  Ref.  8  it  was  pointed  out  that  the  expression  for  the 
viscoelastic  damping  of  extensional  waves  predicts  in¬ 
creased  damping  when  the  extensional  wave  speed  in  the 
constraining  layer  is  less  than  that  in  the  base  plate.  The 
condition  for  optimum  damping.  Appendix  A,  Eq.  (A2) 
appears  in  this  case  as 

g'  =  [(c/c,)^-l]/(l-f-^^)'^^  (5) 

the  denominator  of  Eq.  (A3)  becomes  small,  due  to  the 
term  2  sgn[l  —  (c/cj)^].  This  amounts  to  a  relaxation  type 
of  interaction  between  the  base  plate  and  the  constraining 
layer,  with  strong  attenuation. 
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FIG.  4.  Attenuation  factor  a  relative  to  the  real  pan  of  the  wave  number 
k'  as  a  function  of  frequency  /,  for  a  constrained  niate  in  vacuum,  ac¬ 
cording  to  the  extended  Kerwin  model.  The  physical  parameten  for  the 
three  layers  are  given  in  Table  I.  The  base  plate  is  10-cm  brass,  the 
elastomer  layer  is  1.24  mm  thick.  Solid  curve:  extensional  wave,  con¬ 
straining  layer  2.48-mm  aluminum;  dashed  curve:  extensional  wave,  con¬ 
straining  layer  is  2.48-mm  silver;  dotted-dashed  curve:  flexural  wave,  con¬ 
straining  layer  is  2.48-mm  aluminum. 


As  an  example,  the  aluminum  constraining  layer  is 
replaced  with  a  silver  layer  of  the  same  thickness.  It  is  not 
claimed  that  silver  is  a  good  constraining-layer  material  in 
practice,  but  it  is  used  to  provide  a  demonstration  by 
means  of  an  existing  material.  Silver  was  chosen  because  it 
has  an  extensional  wave  speed  (2933  m/s)  relative  to  that 
of  brass  (3765  m/s)  that  is  about  the  inverse  of  the  exten- 
sional  wave  speed  of  aluminum  (5462  m/s),  relative  to 
that  of  brass.  In  Fig.  4  a  comparison  is  given  for  three 
cases,  based  on  the  extended  Kerwin  model.  In  all  cases  the 
thickness  of  the  brass  base  plate  is  10  cm.  The  solid  line 
gives  the  attenuation  when  the  constraining  layer  is  2.48- 
mm-thick  aluminum.  This  may  be  compared  with  the  at¬ 
tenuation  for  flexural  waves  for  the  same  system,  the 
dotted-dashed  curve.  When  the  constraining  layer  is  2.48- 
mm-thick  silver,  the  maximum  relative  attenuation  in¬ 
creases  by  a  factor  of  about  7,  or  17  dB.  This  is  due  to  the 
ratio  of  the  maximum  of  the  third  factor  in  Eq.  (Al), 
given  in  Eq.  (A3),  for  the  two  cases:  sgn[l  —  (0/03)^]  is 
minus  and  plus,  respectively.  The  second  factor  in  Eq. 
(Al)  is  practically  the  same  in  the  two  cases.  The  sharp 
dip  in  the  dashed  curve  is  a  consequence  of  the  variation  of 
the  extensional  wave  speed  of  the  base  plate  with  fre¬ 
quency,  in  thick-plate  or  exact  theory.  Considering  the 
third  factor  of  Eq.  (Al)  as  a  function  of  c/cy,  it  has  a 
maximum  for  c/Ci  =  [l+g'(\+0^)]‘''^  and  a  minimum 
(zero)  for  c/c3=l.  These  extrema  are  closely  spaced  in 
frequency  for  small  g',  as  is  the  case  in  the  present  example 
(see  also  Appendix  A). 

Figure  5  shows  the  comparison  for  the  aluminum- 
silver  combination  between  the  extended  Kerwin,  hybrid, 
and  exact  models.  One  sees  that  up  to  about  10  kHz  the 
agreement  is  quite  good.  Above  that  frequency,  the  exact 
model  shows  a  peculiar  structure  that  has  not  been  ex¬ 
plained  physically.  The  dip  apparent  in  the  Kerwin  curve  is 
masked  in  the  hybrid  curve  by  the  wider  equivoluminal 
dip,  which  occurs  at  about  the  same  frequency. 
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FIG.  5.  Attenuation  factor  a  relative  to  the  real  part  of  the  wave  number 
fc'  as  a  function  of  frequency  /,  for  a  constrained  plate  in  vacuum.  The 
physical  parameters  for  the  three  layers  are  given  in  Table  I.  The  base 
plate  is  10-cm  brass,  the  elastomer  layer  is  1.24  mm  thick,  the  constrain¬ 
ing  layer  is  2.48-mm  silver.  Solid  curve:  exact  elasticity  theory,  dashed 
curve:  hybrid  model,  dotted-dashed  curve:  extended  Kerwin  model. 

IV.  MAXIMUM  DAMPING  AT  1000  HZ  (CASE  B) 

A.  Constrained  plate  in  vacuum 

In  the  following  examples,  the  thickness  of  constrain¬ 
ing  and  elastomer  layers  are  adjusted  such  that  the  maxi¬ 
mum  viscoelastic  damping  for  extensional  waves  is  near 
1000  Hz.  In  Fig.  6,  the  solid  curve  shows  the  result  for  the 
relative  attenuation  when  the  constraining  plate  is  alumi¬ 
num,  according  to  the  Kerwin  model.  As  before,  for  case 
A,  if  one  replaces  the  aluminum  cover  plate  with  a  silver 
one,  a  higher  attenuation  is  reached.  The  dashed  curve 
shows  the  result  from  the  extended  Kerwin  model,  for  a 
silver  constraining  layer.  The  hybrid  model  gives  the 
dotted-dashed  curve  in  Fig.  6  for  this  case,  which  is  slightly 
different  from  the  extended  Kerwin  model  result. 

The  corresponding  curve  computed  by  the  exact  the¬ 
ory  displays  quite  a  different  picture  above  4000  Hz,  as  the 
solid  curve  in  Fig.  7  shows,  in  comparison  with  the  dashed 
curve  for  the  hybrid  model.  The  first  large  peak  near  5  kHz 


FIG.  6.  Attenuation  factor  a  relative  to  the  real  part  of  the  wave  number 
k'  as  a  function  of  frequency  /,  for  a  constrained  plate  in  vacuum.  The 
physical  parameters  for  the  three  layers  are  given  in  Table  I.  The  base 
plate  is  10-cm  brass,  the  elastomer  layer  is  7.75  mm  thick.  Solid  curve: 
extetided  Kerwin  model,  constraining  layer  13.S-mm  aluminum;  dashed 
curve:  hybrid  model,  constraining  layer  IS.S-mm  silver;  dotted-dashed 
curve:  extended  Kerwin  model,  constraining  layer  15.5-mm  silver. 
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FIG.  7.  Attenuation  factor  a  relative  to  the  real  part  of  the  wave  number 
k'  as  a  function  of  frequency  /.  for  a  constrained  plate  in  vacuum.  The 
physical  parameters  for  the  three  layers  are  given  in  Table  1.  The  base 
plate  is  10-cm  brass,  the  elastomer  layer  is  7.75  mm  thick,  the  constrain¬ 
ing  layer  is  15.5-mm  silver.  Solid  curve:  exact  model;  dashed  curve:  hy¬ 
brid  model. 

corresponds  to  the  first  narrow  peak  in  Fig.  5.  The  char¬ 
acter  of  the  structure  above  10  kHz  may  be  made  compre¬ 
hensible  by  inspection  of  Fig.  8,  where  the  wave  speed  is 
given  as  a  function  of  frequency,  together  with  the  exten¬ 
sional  wave  speeds  for  a  10-cm  brass  plate  and  a  12.6-mm 
silver  plate.  The  wave  speed  in  the  constrained  plate  ap¬ 
pears  to  mostly  follow  the  bra.'ts  curve  up  to  the  point  of 
intersection  of  brass  and  silver  curves.  From  there  it  fol¬ 
lows  the  silver  wave  speed  for  a  while,  and  then  meanders 
in  the  space  between  the  silver  and  brass  curves,  to  rejoin 
the  brass  curve  at  the  next  point  of  intersection.  Fragments 
of  other  branches  (not  shown)  were  found,  but  were  not 
followed  throughout  the  frequency  range. 

B.  Fluid-loaded  constrained  plate 

Figure  9  depicts  the  results  according  to  the  hybrid 
model  for  a  constrained  plate  loaded  by  water  on  both 
sides,  on  the  side  of  the  constraining  layer,  and  on  the 
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FIG.  8.  Phase  speed  c  as  a  function  of  frequency  /  for  a  constrained  plate 
in  vacuum,  compared  with  phase  speeds  for  single  plates,  according  to 
exact  elasticity  theory.  The  physical  parameters  for  the  three  layers  are 
given  in  Table  I.  Solid  curve:  base  plate  is  10-cm  brass,  the  elastomer  layer 
is  7.75  mm  thick,  the  constraining  layer  is  15.5-mm  silver;  dashed  curve 
#1:  10-cm  thick  brass  plate  in  vacuum;  dashed  curve  #2:  15.5-mm-thick 
silver  plate  in  vacuum. 
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FIG.  9.  Attenuation  factor  a  relative  to  the  real  part  of  the  wave  number 
A:'  as  a  function  of  frequency  /,  for  a  constrained  plate,  according  to  the 
hybrid  model.  The  physical  parameters  for  the  three  layers  are  given  in 
Table  I.  The  base  plate  is  10-cm  brass,  the  elastomer  layer  is  7.75  mm 
thick,  the  constraining  layer  is  IS.S-mm  silver.  Solid  curve:  water  on  both 
sides;  short-  and  long-dashed  curve:  water  on  side  of  constraining  layer  or 
opposite  side  (these  cases  are  indistinguishable  at  this  scale);  dashed 
curve:  no  fluid  loading. 

opposite  side.  The  latter  two  cases  are  indistinguishable  at 
the  scale  of  the  figure.  The  dashed  curve  is  for  an  unloaded 
composite  plate.  One  may  clearly  distinguish  a  lower  fre¬ 
quency  range,  where  viscoelastic  damping  dominates,  and 
an  upper  range,  where  radiation  is  the  main  damping 
mechanism,  and  where  the  curves  look  very  similar  to 
those  of  Fig.  1  for  a  single  plate.  The  exact  model  for  fluid 
loading  on  both  faces  shows  a  slightly  different  picture 
(Fig.  10):  a  third  peak  appears  between  the  two  peaks  of 
Fig.  9,  obviously  related  to  the  high  peak  identified  in  Fig. 
7  near  5  kHz,  and  thus  indicative  of  viscoelastic  damping. 
It  is  unexpected,  though,  that  the  radiation  peak  near  20 
kHz  is  lower  than  predicted  by  the  hybrid  model  (dashed 
curve),  which  coincides  at  high  frequency  with  the  radia¬ 
tion  damping  by  a  single  plate  (dotted-dashed  curve).  This 
would  seem  to  be  another  manifestation  of  interaction  of 
radiation  and  viscoelastic  damping. 


FIO.  10.  Attenuation  factor  a  relative  to  the  real  part  of  the  wave  number 
k'  as  a  function  of  frequency  /,  plate  loaded  by  water  on  both  sides.  The 
physical  parameters  for  the  three  layers  are  given  in  Table  I.  The  base 
plate  is  lO-cm  brass.  Solid  curve  (exact  model)  and  dashed  curve  (hybrid 
model)  for  constrained  layer,  with  7.73-mm-thick  elastomer  layer  and 
15.5-mm  silver  constraining  layer;  dotted-dashed  curve  (exact  model)  for 
single  plate. 


mrooENCV  (mi) 

FIG.  11.  Attenuation  factor  a  relative  to  the  real  part  of  the  wave  number 
A'  as  g  function  of  frequency  /,  for  a  constrained  plate,  according  to  the 
exact  model.  The  physical  parameters  for  the  three  layers  are  given  tn 
Table  I.  The  base  plate  is  lO-cm  brass,  the  elastomer  layer  is  7.75  mm 
thick,  the  constraining  layer  is  15.5-mm  silver.  Solid  curve;  both  sides 
loaded  by  water;  dashed  curve:  side  of  constraining  layer  loaded  by  water; 
dotted-dashed  curve:  opposite  side  loaded  by  water. 

In  Fig.  1 1  the  effect  of  fluid  loading  is  demonstrated, 
for  fluid  loading  on  both  sides,  on  the  side  of  the  constrain¬ 
ing  layer,  and  on  the  opposite  side,  as  computed  by  the 
exact  model.  The  result  for  the  fluid  loading  on  the  side  of 
the  constraining  layer  is  strikingly  different  from  the  two 
other  fluid-loading  modes;  it  appears  to  reflect  the  behavior 
for  the  unloaded  plate  of  Fig.  7,  with  an  increasing  radia¬ 
tion  contribution  added  for  increasing  frequency.  One  is 
inclined  to  interpret  this  as  an  effect  of  the  motion  in  the 
plate  being  drawn,  at  high  frequencies,  to  the  side  of  the 
fluid  loading  the  plate;  if  the  fluid  is  at  the  constraining- 
layer  side,  both  viscoelastic  and  radiation  damping  effects 
are  visible.  Thus,  the  situation  is,  in  a  sense,  complemen¬ 
tary  to  that  for  flexural  waves,  but  caution  is  required, 
since  the  structure  of  the  curves  is  more  complicated  for 
extensional  waves. 


V.  SUMMARY  AND  CONCLUSIONS 

A  thin-plate  theory  for  extensional  waves  in  a  fluid- 
loaded  plate  was  presented  that  gives  good  agreement  with 
the  exact  theory  in  the  case  of  loading  by  water  on  both 
sides. 

The  damping  of  extensional  waves  in  infinite  con¬ 
strained  plates  without  and  with  fluid  loading  was  analyzed 
by  three  models:  the  extended  Kerwin  model,  the  hybrid 
model,  and  the  exact  model. 

The  extended  Kerwin  model  is  of  importance  for  de¬ 
signing  constrained-layer  damping  applications  with  spe¬ 
cific  properties,  since  it  expresses  the  attenuation  constant 
explicitly  in  terms  of  physical  and  geometric  parametei'S.  It 
is  limited  to  cases  where  the  elastomer  and  constraining 
layer  are  thin  relative  to  the  base  plate,  and  it  docs  not 
account  for  fluid  loading. 

The  hybrid  model  is  useful  for  predicting  the  effect  of 
fluid  loading,  and  is  reliable  for  most  frequency  ranges  of 
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practical  imptortance.  It  requires  fewer  operations  than  the 
exact  model.  Both  the  hybrid  and  the  exact  model  imply 
finding  the  complex  roots  of  an  equation  formed  by  setting 
the  determinant  value  of  a  matrix  equal  to  zero.  For  the 
hybrid  model,  this  matrix  has  order  4x4,  5x5,  or  6X6, 
depending  on  whether  the  constrained  plate  is  in  vacuum, 
loaded  by  fluid  on  one  side,  or  two  sides,  respectively.  For 
the  e.xact  model  the  order  is  12X  12,  13x  13,  or  14x  14. 

The  exact  model  is  the  best  model  within  the  limita¬ 
tions  of  linear  structural  acoustics.  It  is  recommdended  as 
a  backup  to  the  two  other  models,  especially  in  situations 
that  differ  from  the  traditional  CLD  conditions  of  metallic 
plates  with  thin  added  layers  at  modest  frequency- 
thickness  products. 

The  constrained  layer  damping  of  extensional  waves 
displays  a  few  features  not  applicable  to  the  damping  of 
flexural  waves. 

(i)  It  is  necessary  to  use  a  constraining  layer  of  a 
material  different  froni  the  base  plate  material. 

(ii)  There  is  a  peculiar  relaxation-type  behavior  of  the 
damping  when  the  extensional  wave  speed  in  the  constrain¬ 
ing  layer  is  less  than  that  in  the  base  plate. 

(iii)  For  one-sided  fluid  loading,  there  are  indications 
that  the  wave  amplitude  on  the  loaded  side  increases  with 
increasing  frequency,  contrary  to  the  case  of  flexural 
waves,  where  the  wave  energy  moves  to  the  face  opposite 
to  the  loaded  face. 
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APPENDIX  A:  EXTENDED  KERWIN  MODEL 

In  Ref.  2  the  following  expression  was  proposed  as  an 
extension  of  the  Kerwin  model: 


For  high  frequency,  the  extensional  wave  speed  c  in 
the  base  plate  decreases  rapidly  as  a  function  of  frequency. 
For  the  case  where  the  thin-plate  extensional  wave  speed  of 
the  constraining  layer  is  less  than  that  in  the  base  plate,  the 
third  factor  of  Eq.  (A1 ),  considered  as  a  function  of  c/c,, 
has  a  minimum  for  c  =  Cj,  equal  to  zero,  and  a  maximum 
for  (c/Ci)'=\-tg'{l  equal  to  g'(  1 +^^)/(4-f 

The  minimum  and  maximum  are  very  close  together,  since 
g'  is  small  in  this  frequency  range. 


APPENDIX  B:  MODEL  FOR  EXTENSIONAL  WAVES  IN 
A  FLUIJ-LOADED  THIN  PLATE 

Thin-plate  theory  for  flexural  waves  in  a  fluid-loaded 
plate  may  be  found  in  Ref.  5.  In  that  reference  a  term  q  is 
introduced  into  the  equation  for  a  beam'^  as  a  force  per  unit 
length  and  by  analogy  appears  in  the  equation  for  a  plate  as 
a  force  per  unit  area,  equal  to  the  pressure  of  the  fluid.  This 
leads  to  the  following  wave  equation  for  flexural  waves  <n  a 
fluid-loaded  plate^ 

.  d^w 

DV*w+ psh -^= .  (Bl) 

where  V*  is  the  biharmonic  operator,  D  is  the  bending 
stiffness  of  the  plate,  D—{\/\2)Eh^/(  \  —  v^),  Fis  Young’s 
modulus,  h  is  the  plate  thickness,  v  is  Poisson’s  ratio,  w  is 
the  displacement  perpendicular  to  the  plate,  p,  the  density 
of  the  plate  material,  and  p  the  pressure  in  the  fluid  (fluid 
on  one  side).  For  two-sided  fluid  loading  the  pressure  term 
is  doubled. 

The  goal  of  the  present  derivation  is  to  arrive  at  a 
similar  equation  for  extensional  waves  in  a  thin,  fluid- 
loaded  plate,  in  the  form 

E  d^w 


a  _  g'[l-(c/c,)^]^ 

E\h,  [  1  -  (c/c3)'+g'l4^'g''  ’ 


(Al) 


where  a/k’  —  ihe  attenuation  constant  a  relative  to  the  real 
part  of  the  wave  number  k'\  ^=the  loss  tangent  of  the 
elastic  shear  modulus;  c  =  the  speed  of  flexural  or  exten¬ 
sional  waves  from  the  thick-plate  or  exact  theories;  Cj^the 
extensional  wave  speed  in  layer  3;  F]  3  =  the  Young’s  mod¬ 
ulus  of  layers  1  and  3;  F^=F/(1— v^),  where  v  is  the 
Poisson  ratio;  g=the  shear  parameter,  g—Gi/ik'^E^hjh^), 
and  g=g'(  1  +</3):  <7j  =  the  complex  shear  modulus  of  the 
elastomer,  Gi  =  (/jf  1  -f  iP)  \ and  A,  — the  thickness  of  layer 
i,  t=  1,2,3;  r—l/A  for  flexural  waves  and  1/2  for  exten¬ 
sional  waves. 

The  third  factor,  considered  as  a  function  of  the  real 
part  of  the  shear  parameter  g',  has  a  maximum  for 

g'=|l-(c/c3)2|/{l+/3^)''^  (A2) 

the  maximum  value  of  this  factor  is 


11-(C/C3)'1 

2(  I  -f/3^ ) 2  sgn ( 1  -  (c/cIFI  ’ 


(A3) 


where  F  is  a  fluid  loading  term,  with  the  dimension  of  a 
force  per  unit  length  cubed.  There  was  no  indication  found 
in  the  literature  how  this  term  could  be  related  to  the  fluid 
pressure  from  first  principles  in  the  case  of  extensional 
waves,  in  contrast  to  flexural  waves  as  shown  in  Eq.  (Bl ). 

Therefore,  the  problem  was  approached  by  means  of 
thick-plate  theory  for  extensional  waves.'*  Thick  plate 
theory  carries  the  approximation  implied  by  thin-plate  the¬ 
ory  to  the  next  power  of  the  dimensionless  wave  number 
kd,  where  k  is  the  wave  number  and  d  the  half-thickness  of 
the  plate.  The  result  for  flexural  waves  is  known  as 
Timoshenko-Mindlin  plate  theory."  In  this  type  of  theory 
the  fluid  loading  appears  naturally  as  a  consequence  of  the 
introduction  of  the  structural  equations  for  the  plate.  Once 
the  fluid-loading  term  has  been  introduced,  thin-plate  the¬ 
ory  with  fluid  loading  follows  by  taking  the  limit  kd  -0. 

By  integrating  the  differential  equations  of  motion 
from  exact  elasticity  theory  or  the  corresponding  momen¬ 
tum  equations  in  the  direction  z  perpendicular  to  the  di¬ 
rection  of  propagation  x  (parallel  to  the  faces  of  the  plate), 
one  obtains  the  following  structural  equations  for  exten¬ 
sional  (symmetric)  waves:'*’ 
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c^U  dy 

(A  +  2G)  ^+/l™  +  pyC/  =  0  (83) 

and 

2  ,  ,  dU  2  ,  , 

-  ^/VG ^-2{A  +  2G);^^-2A  — +  ^^+3  rf'A'-O, 

(84) 

where  the  displacement  u  in  the  direction  of  propagation  x 
is  given  by  a  z-independent  function  U(x),  while  the  z 
displacement  w  is  approximated  by  an  antisymmetric  lin¬ 
ear  expression  in  z,  w—zy(x).  Here  A  and  G  are  the  Lame 
constants,  o)  is  the  angular  frequency,  d(=h/2)  is  the  half¬ 


thickness  of  the  plate,  and  k  is  a  coefficient  accounting  for 
the  varying  shear  in  the  plate.  The  fluid  loading  is  repre¬ 
sented  by  the  term  the  superscript  .v  indicates  exten- 
sional  waves,  for  which  (f  —  a^  +a  ,  where  are  the 
normal  stresses  at  the  two  faces  of  the  plate,  equal  to  the 
negative  of  the  pressures  on  the  faces. 

Assuming  a  harmonic  straight -crested  wave  with 
space  and  time  dependence  according  to  exp /(«/  — Ax), 
where  k  is  the  (complex)  wave  number,  k~k'—ia,  one 
finds  the  dispersion  relation  by  equating  to  zero  the  deter¬ 
minant  of  the  matrix  of  the  coefficients  of  the  ampli¬ 
tudes  U,  y  in  the  two  Eqs.  (83)  and  (84) 


-A^(A+2G)-|-py  -lAA 

2/AA  -  IcP-k^K^G— 2  ( A  -f  2G)  -f  +  q'/y 

By  manipulating  the  equation  Det(Mj)  =0,  one  finds 
,  (A  +  2G)~p.,c^  A 

A  {k  +  lG)+\{kd)\^G -\p^{kd)^-(f /(2y) 


=  0. 


(85) 


(86) 


The  thin-plate  model  follows  by  dropping  terms  of 
order  {kd)^.  By  expanding  the  determinant  one  obtains 

(A-f2G)^-p^^(A-f2G)-A^ 

-[{MG)-pf^W/{2x)=0.  (87) 

If  one  divides  this  equation  by  (A-I-2G)  =pjC^  and  applies 
(A-I-2G)— AV(A-i-2G)=£/(1  — v^),  there  results  the 
equation 

£/(  1  -  v^)  -p/-  [( 1  - {c/Cj)^W/{2x)=^0,  (88) 

where  cj  is  the  dilatational  wave  speed  in  the  plate  mate¬ 
rial. 

To  see  how  the  fluid  loading  appears  in  the  wave  equa¬ 
tion  for  w(—xd),  the  z  displacement  of  the  plate  surface  at 
the  positive  side,  one  compares  this  expression  with  the 
proposed  thin-plate  equation  with  fluid  loading  Eq.  (82) 
Comparing  Eqs.  (82)  and  (88),  one  sees  that  the  fluid¬ 
loading  term  F  is  given  by 

F=l\-ic/cj)^](kd)^q‘/(2d).  (89) 

The  form  of  this  equation  suggests  that  it  would  be 
difficult  to  derive  this  expression  from  mechanical  princi¬ 
ples. 

To  obtain  an  equation  in  the  variable  k  one  uses  Eul¬ 
er’s  equation.  For  fluid  loading  on  both  sides,  one  has  (f 
=  2&»^Pow/(A:^  — where  po  is  the  density  of  the  me¬ 
dium  and  Aq  the  wave  number  in  the  medium.  This 
expression  is  inserted  into  Eq.  (88).  Similar  to  the  devel¬ 
opment  in  Ref.  5,  the  resulting  equation  acquires  a  simple 
form  if  one  introduces  a  new  variable  r  defined  by  (Aqt)^ 
=  A^  — Aq.  After  some  manipulation,  the  result  is  a  third- 
degree  equation  in  r 
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T^-a7^  +  r[\-{co/Cp)^]~n[l-(co/Cj)-]=0, 

(89) 

where  Cp  is  the  wave  speed  for  extension?!  waves  in  a  thin 
plate.  A  characte,isac  frequency  Uj.  and  a  dimensionless 
frequency  fi  are  introduced  by  ci>c=ip/po)cp'(C(pl)  and 

il  =  U)/ £t)j.. 
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